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ABSTRACT 

A pole-count analysis of the infrared 2MASS survey is presented, in order to identify 
faint stream-like structures within the halo of the Milky Way. Selecting stars with 
colours consistent with M-giant stars, we find a strong over-density of sources on a 
stream with pole {i = 95°, 6 = 13°), which corresponds to the pole of the orbit of 
the Sagittarius dwarf galaxy. This great-circle feature of width ^ 12°, contains ^ 5% 
of the late M-giants in the Halo. No other stream-like structures are detected in M- 
giants in the 2MASS Second Incremental Data Release (2IDR), and in particular, 
we find no evidence for a stellar component to the Magellanic Stream. This suggests 
that the present accretion rate of low-mass satellites with a luminous component 
is very low, and the formation of the luminous component of the Halo must have 
been essentially complete before the accretion of the Sagittarius dwarf galaxy, more 
than 3 Gyr ago. We also search for great-circle streams using almost all high-latitude 
> 30°) sources in the 2IDR dataset. No narrow great-circle streams of width 
0.5°-2° were found, though we were only sensitive to relatively nearby (< 17kpc) 
remnants of massive (10^ Mq) globular clusters. If the Galactic potential is close to 
being spherical, as some recent observations suggest, the lack of observed great-circle 
streams is consistent with the presence of dark matter substructures in the Halo. 
Although alternative explanations cannot be ruled out from our analysis of the 2IDR 
dataset, future experiments with better statistics have the potential to reveal the 
heating effect of dark matter substructure on stellar streams. 

Key words: halo — Galaxy: structure — Galaxy 



1 INTRODUCTION 



Accord 



iig Co out cmreuC uiideisCaiidiiig of Che fomiaCiuii of 



galaxies, galaxy hales build up over time by the conti nued 
accreti(in o f smallm galacLic ui sub-galacLic sLmcCmus (6rfe~ 



et al. 1S?94 ). Of the theori e s that quantify this formation pro- 



cess, the Cold Dark Matter theory gives at present the best 
explanation of available data, but there are several inconsis- 
tencies ilii Lilly plcLuiy LhaC leiuaiu Lu by explained. Possibly 
the mo st troublesome of these is the inconsistency between 
the predicted and observed density profiles of galaxy ha- 
los. CDM predicts strongly centrally-peaked (p oc r^^ to 



* This I lublication makes use of data products from the Two Mi- 
cron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



p oc r ^•^) mass distrib utions (Navarro, Frenk & White 
1997; Moore et al. 1999a), contrary to what is inferred from 



the slowly-rising rotation cu rves of highly dark matter dom- 



inated dwarf disk galaxies (^ote et al. 2000 



ie Blok et al 



2001; Marchesini et al. 2001). Another problem is the ap- 
parent disagreement between the expected number of sub- 
structures in CDM halos compar ed to the visible substruc- 
tures in the halos o f real galaxies ( Klypin et al. 1999| ; 



Moore 



et al. 1999b, 2001). It appears that a proper treatment of 
the ionizing effect of the intergalactic UV background, which 
disperses the baryonic component of low mass proto-galaxies 
into the intergalactic medium, can explain the absence of co- 



pious low-mass substructures in galaxy halos ( Bullock et al 



200C| ; pomerville 200l| [Tully et al. 20011 ). There is also some 
concern that current cosmological simulations are of insuffi- 
cient resolution: for instance, one needs excellent numerical 
resolution to properly resolve the wakes in the dark mat- 
ter distribution which give rise to dynamical friction. If the 
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resolution is not high enough, dynamical friction will be ar- 
tificially low, and stru ctures will take lon ger to merge than 
they would in reality (van Kampen 200C). 

Further observations are required to help resolve these 
issues with CDM or to guide galaxy formation theory. Data 
that constrain the structure of the luminous and dark mat- 
ter in the inner regions of galaxies are relatively easy to 
obtain, but in those regions the dynamical times are short 
and much of the detailed information about the formation of 
the galaxy has been phased-mixed and is now largely lost. 
In the outer regions of galaxies, typical dynamical mixing 
timescales are long, and so valuable information on the for- 
mation of the galaxy can be recovered, if only we can find 
tracer material or dynamical probes. This is, of course, the 
challenge of studying the halos of galaxies: their very low 
surface brightness renders them unobservable, in general, to 
present instrumentation. 

One solution to this problem is to study the halo of the 
Milky Way and other nearby galaxies, where the halo stars 
are resolved. Radial velocities can be obtained in this situ- 
ation to complement the starcounts data, allowing detailed 
dynamical modeling of the stellar population. In this way 
one cal l studv the structure and substructure of the under- 



lying d< rk matter even if it is not traced by stars. The nearby 
universe therefore gives us a unique and more complete view 
of the dark matter, which we cannot gain from statistical 
studies of more distant systems. The richness of the infor- 
mation in this region will be even more striking after the 
completion of the astrometric satellite missions GAIA and 
SIM: we will then also have precise measurements of the 
transverse motion of a huge sample of Galactic and Local 
Group stars. 



In a companion paper (Ibata, Lewis & Irwin 2002 
Paper 1), we discuss the possibility of using narrow tidal 
streams from defunct globular clusters to ascertain the 
lumpiness of galactic halos. We show that the heating due 
to repeated close encounters with the dense dark matter 
clumps predicted by CDM and other structure formation 
theories spread an otherwise narrow tidal stream into a 
coarse band. Individual streams from the expected popu- 
lation of tidally destroyed globular clusters should be read- 
ily detectable with the GAIA all-sky survey and the effect 
of the broadening of the globular-cluster tidal stream should 
be a strong effect if CDM predictions are correct. The GAIA 
dataset will be by far the most efficient way to reveal these 
ancient tidal streams, since there will be proper motion mea- 
surements, radial velocities as well as a good spectral energy 
distribution from multi-band photometry. The kinematic 
data will allow selection of stars in a 5-dimensional phase- 
space (perhaps even in 6-dimensions for nearer streams for 
which parallax measurements are possible). 

However, most of the giant-branch stars of the streams 
that will be detected by GAIA have already been cataloged 
in current all-sky surveys. Ancient accretions of low-mass 
galaxies and globular clusters in the outer halo of the Milky 



Way sh ould give rise to long stellar streams (Johnston, Hern- 



quist & Bolte 1996). These streams may be detected by 
searching for long band-like overdensities on the sky if the 
Halo potential is not very flattened. 

In this contribution we undertake a preliminary analysis 
of the 2 Micron All Sky Survey (2MASS) to investigate if 
any streams can be identified from sky position together 



with simple colour-magnitude cuts that reject foreground 
contamination. Our two-dimensional search is therefore only 
likely to recover the structures of greatest contrast over the 
background. 

Some tidal streams from disrupting stellar systems have 
already been identified. Leon, Meylan & Combes (2000) have 



found low-level tidal streams around several globular clus- 
ters, and another very nice example has been reveale d by the 
Sloan Digital Sky Survey (Odenkirchen et al. 2001). These 



globular cluster streams, however, are a most a few degrees 
long and it has not been possible to trace them over a large 
extent on the sky, which is needed if we are to be able to 
probe the lumpiness of the halo dark matter. Large-scale 
tidal streams of more massive systems such as disrupted 
dwarf galaxies have also b een found. In another companion 
paper (Ibata et al. 2001c), we analysed the distribution of 
Carbon stars (C-stars) in the halo of the Milky Way, and 
found that more than half of these are distributed along a 
great circle on the sky. These C-stars trace the tidal stream 
of the Sagittarius dwarf galaxy, which has been disrupted 
over the course of its many close encounters with the Milky 
Way. In extragalactic sys tems, huge tidal st reams hav e been 
found aroun d NGC 5907 (|Shang et al. 1998|) and M31 ([ibata 



et al. 2001b) 



2 SAMPLE 1: SELECTION OF LATE M-GIANT 
HALO STARS 

The dataset we use is the Second Incremental Data Release 
(2IDR) of the 2MASS project. [| This release covers 47% 
of the sky in a rather patchy fashion, as can be seen from 
Figure 1, which shows the density of point-source detections 
with K < 14.3 (the 10 sigma limit of the 2MASS photom- 
etry). The vast majority of the stars in Figure 1 belong to 
the Galactic disk and bulge; these outnumber spheroid stars 
by several orders of magnitude. They are the background 
"noise" that has to be removed, or accurately modelled out, 
from the sample. 

Since the structures we are interested in revealing are 
at large distances in the Halo, it is convenient to remove 
first the disk stars from the 2MASS sample. To do this 
we took advantage of the structure of the J — H , H — K 



colour-colour diagram. As shown by Bessel & Brett (1988) 



dwarf and giant stars beyond Ml separate cleanly, the M- 
giants being 0.1 to 0.3 magnitudes redder in J — H than 
M-dwarfs at a given H ~ K colour. Selecting stars with 
(J -//)(, > 0.75 + 2S{J - H), {H-K)o < 1.5-25{H-K), 
and (J - H)o > {H - K)o + 0.4 + 2S{H - K) removes 99.9% 
of all of the 2IDR 2MASS point-source detections, the ma- 
jority of which are nearby disk dwarfs. (The colours were 
corrected for interstellar reddeni ng by interpolating i n the 
IRAS-DIRBE reddening map of ^chlegel et al. (199^ ), us- 
ing Aj = 0.902£'(B - V), Ah = 0.576E{B - V), Ak = 
0M7E{B - V)). Though the relative fraction of dwarf and 
giant stars in the disk and spheroid may be similar, by se- 
lecting only giant stars, we are able to clean much of the 
disk contribution to the high-latitude sky, but even with 



^ see ^ittp : //www, ipac . caltech.edu/2mass 
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Table 1. Parameter cuts employed for selection of the M-giant 
star sample. 
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< 
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this nlpaning wp find that thp. mn tarn in at ion at ?n-.Sn Hp- 
grees f jom the Plane is still sufficiently high that we are 



restricted to using a sample with |b| > 30°. 

Given a characteristi c absolute magnitude of RG B M- 
giants of Mr ^ -5.5 ( [Nikolaev fc Weinberg 200(]| ), the 
bright limit of isT = 8.0 of the 2MASS data corresponds to a 
lower distance limit of 5 kpc for the giant stars. This distance 
is more than three thick disk scale-heights (zh ~ 800 pc 
Reyle fc Robin 200 1[ ) for fe > 29°. (The 2MASS faint limit 
of K — 14.3 corresponds to an upper distance limit of 
- 100 kpc.) 

We decided to ignore the data at Galactic latitudes 
lower than |fe| = 30°, as contamination from the Galactic 
thick disk begins to be significant below that latitude. Strong 
reddening can move dwarfs in the region of the colour-colour 
diagram inhabited by giants, so we also excised from our 
final dataset all regions of the sky wh ere the reddening de - 
duced from the IRAS-DIRBE map of ^chlegel et al. (199^ ) 
exceeded E{B — V) — 0.1. Because of this reddening cut, 
and to maintain a homogeneously deep map, we further re- 
stricted the K-band limit to K = U.3-0.367E{B-V)max = 
14.26. At Galactic latitudes |&| > 30°, the area of sky which 
is removed by this extinction cut is very small (0.6%) and is 
confined to regions around nearby molecular clouds. Finally, 
two more regions are cut out of the dataset: a circular region 
of radius 10° around the LMC and SMC. The resulting map 
is shown in Figure 2, which shows the 1106 point-sources on 
the sky that remain after these various cuts; the area of sky 
available for analysis (see Figure 3) is 26.4%. 



3 POLE-COUNT ANALYSIS OF SAMPLE 1 

One way to identify long band-like structures is through a 
pole-count analysis. The pole-count analysis is very simple: 
the number of sources situated within a band centered on 
a great circle on the sky is counted up, and assigned to the 
corresponding pole. This is repeated for all possible poles. In 
this way, a band-like structure on the sky becomes two peaks 
in the pole-counts diagram (two peaks due to the symmetry 
of the mapping), whereas a localised peak in counts on the 
sky becomes two bands in the pole counts diagram. Due to 
the complicated irregular nature of the window function, it 
is more convenient to perform the pole-counts on binned 
data: to this end we binned up the sky into a rectangular 
3600 X 1800 grid. To correct for the sky regions for which 
we have no data, we add in the average number density of 
sources into each pixel in the image (see Figure 3). 

The pole-counts analysis is particularly powerful if the 



underlying potential in the region of the Galaxy inhabited 
by the stream is spherical, in that case the stellar stream 
does not precess and it remains on a great circle on the sky 
forever. In a halo with a flattened potential, differential pre- 
cession will eventually phase-mix the stream, thereby filling 
the box-orbit. In this situation, only strong recent accretions 
would be detected by the pole-counts technique. 

Fortunately for the pole-counts method, the outer re- 
gions of the Milky Way appear to have an almost spheri- 
cal potential. The evidence for this comes partly from the 
great-circle distribution of C-s tars from the tida. 1 stream of 
the Sagittarius dwarf galaxy (Ibata et al. 2001c), and also 
from the distribution of metal-poor spheroid star s towards 
R ^ 20 kpc of the inner Halo region surveyed by Chiba & 



Beers (2001) 



The free parameter that we need to choose is the width 
of the equatorial bands. Clearly, narrow bands are better 
suited to identify narrow structures. We chose not to probe 
band widths below 1°; this was motivated by the simulations 
of Paper 1, where model globular clusters gave rise to stream 
widths of ^ 2 kpc (equivalent to ~ 1° at 100 kpc). The anal- 
ysis was repeated with band widths 1°, 1.5°, 2°, . . . , 20°, and 
peaks > 3(T above the background were identified in each 
polecounts image. The result of the pole-count for a band 
width of 12.0° (which gives the highest S/N peak) is shown 
in Figure 4. Of our Halo sample of 1106 sources, 15% are lo- 
cated within this 12° band, compared to 10% that would be 
expected from a Poisson distribution. How significant is this 
peak given the strange window function of the 2IDR data? 
To answer this question, we have performed 1000 Monte- 
Carlo simulations, drawing in each case 1106 random points 
over the sky within the data regions displayed in Figure 3. 
The analysis was repeated in exactly the same way as for 
the real data, searching for peaks in the pole-counts dia- 
gram. In the 1000 simulated datasets, the highest deviation 
was a 12.9% peak, found at an {£, b) location far from the 
Sgr pole. Interpreting this maximum random deviation as a 
3.2(7 peak, implies that the peak in the pole-counts distri- 
bution detected in Figure 4 has a significance of ~ 5.5cr and 
is therefore clearly real. 



4 SAMPLE 2: ALMOST NO RESTRICTIONS 

In selecting only late M-giant stars from the 2MASS dataset, 
we previously retained only a small fraction of the full sam- 
ple of Halo stars observed by that survey. To try to overcome 
this handicap, we undertook a second analysis, this time 
with a much broader colour selection, aimed at including K- 
stars and early M-stars as well as late M-stars selected for 
the previous sample. The selection criteria are simple colour 
limits: -0.1 < [H - K)o < 0.5, 0.1 < (J - H)o < 1.0; a 
faint magnitude limit: Kq < 14.26; and a selection of only 
low extinction sky E{B — V)< 0.1. Selecting from the high 
latitude sky (|6| > 30°) in the 2IDR with these criteria gives 
a much larger sample of 5.9 x lO'' point sources. The chal- 
lenge now is to be able to extract a stream population with 
at most a few thousand sources from such a numerous fore- 
ground population. 

To successfully identify streams with the polecounts 
method, we need to model out these contaminants. How- 
ever, the distribution of Galactic sources is not uniform, and 
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it is only approximately symmetric about the Galactic ma- 
jor and minor axes. We found that it was necessary to create 
a model of the smooth Galactic foreground which was able 
to follow the variations on a scale of ~ 10° ; smoother models 
give rise to unacceptably large residuals which dominate the 
subsequent polecounts. The complex 2IDR window function 
also makes it difficult to fit an accurate model: e.g., poly- 
nomial fits of high enough order to follow the ~ 10° Galac- 
tic variations tend to blow up in the unconstrained regions, 
leaving systematic residuals around the holes. 

By trial and error, we managed to construct an accept- 
able model of the Galaxy by smoothing the data. Ifowever, 
to avoid edge effects around the holes in the 2IDR dataset, 
the holes needed to be filled before this smoothing opera- 
tion. The empty regions were filled preferentially with data 
from the symmetric pixel refiected about I = 0, otherwise 
by data from the pixel refiected about 6 = or the pixel 
refiected about I = and b = 0. After this operation, 8.6% 
of the sky still lacked data; these remaining holes were filled 
with a Poisson random deviate drawn from a third order 
polynomial fit to each longitude line separately. The result- 
ing all-sky image was smoothed with a 5.1° bi- linear kernel 
(taking the average of counts within ±4(7 of the median). 

An image of the difference between the original dataset 
and the smoothed model was created and the window func- 
tion was reinserted into the image. This difference image is 
now flat to ~ 2%. The difference image was smoothed again 
with a 5.1° bi- linear kernel, but this time, those regions with 
data on either side of a hole were (for the purposes of the 
smoothing) made contiguous. After removing this smoothed 
model, the final dataset is flat to ~ 0.5%. 

This image of the residuals of the 2IDR from a 
smoothed model is appropriate for searching for narrow 
great-circle streams of width up to ~ 2°. Polecounts di- 
agrams were computed for bands of 2.0°, 1.0° and 0.5°. 
Within a 2° wide band (1.7% of the sky) 10^ stars are ex- 
pected, on average, with la variations of 320 counts. The 
corresponding polecounts diagram had an R.M.S. dispersion 
450 counts, close to this Poisson limit. Since 2° wide bands 
will give rise to ~ 2° wide peaks in the polecounts image, 
there are ~ 5000 "meta-pixels" in that image, so > 4a peaks 
are required to reject most false detections due to random 
noise. There is a single peak with more than 1800 counts, 
which appears as an arc in the polecounts diagram near 
{£ = 0°,b = 0°). The feature is clearly not real, however, 
as the corresponding equatorial band passes through a long 
narrow gap in the window function that skims the outer 
bulge of the Milky Way (near b — —30°). The polecounts 
peak disappears if we remove this small region of sky from 
the analysis, so the peak is clearly an artifact of insufflcient 
accura cy in the smooth Galaxy model in that region. 



pol e-counts diagram o f Figure 4 is very similar to that found 
by Ibata et al. (2001c) from Halo C-stars. M-giants are po- 



With very narrow 0.5° bands, the R.M.S. dispersion of 



the polecounts image was 215 counts, so streams with a few 
as 860 stars could, in principle, have been detected. No such 
peaks were detected, however. 



5 DISCUSSION 

5.1 Wide streams from disrupted dwarf galaxies 

The 2MASS 2IDR data reveal a strong overdensity along 
the great circle track of the Sgr dwarf galaxy. Indeed, the 



tentially older stars than C-stars, and as such they can be 
useful in the context of analysing the Sgr stream, allowing 
us to probe the stars which were stripped from the dwarf 
galaxy at earlier stages in the disruption. To constrain this 
possibility, the radial velocities of the M-giants needed to be 
measured. 

For bands of width < 2°, the Sgr polecounts peak breaks 
up into several lower S/N ~ 3 peaks, presumably because 
the orbit is not precisely on a great circle. This could poten- 
tially be due to a fiattening of the Galactic potential (the 
Galactic disk clearly provides non-radial forces); an analy- 
sis of this effect could be attempted with the full 2MASS 
dataset. 

It is surprising not to have found stream-like structures 
in the Halo from other dwarf satellites, such as the Mag- 
ellanic Clouds. Our previous polecount analysis of Halo C- 
stars provided a hint that the Magellanic Stream had been 
detected. Despite the better statistics, we have not been 
able to confirm that tentative detection. It is possible that 
this refiects a difference in stellar populations between the 
Sagittarius and Magellanic streams, the latter possibly being 
richer in C-stars due to younger age. Another explanation is 
that the current sky coverage excludes most of the region in 
and around the Magellanic Clouds, this makes the detection 
of the Magellanic Stream unlikely with the current dataset 
if the stars are distributed in a similar fashion to the gas in 
the Magellanic Stream. 

Our models of the disruption of the Sgr dwarf suggest 
that the streams are very long lived, remaining coherent 
for a large fraction of a Hubble time. However, differential 
precession can cause a stream to smear out over the sky, until 
the box orbit is filled up entirely. Two possibilities present 
themselves. Either the Halo is fairly flattened, contrary t o 
the results of[lbataet al. (2001a|) and|Chiba fc Beers (200]]). 



In this case streams precess quickly, and pole-counts give us 
a snapshot of accretions that have taken place only very 
recently, on a timescale of ~ 1 Gyr or less (depending on the 
orbit). If the Halo is close to being spherical, polecounts give 
us access to much more ancient accretion events. Without 
a conclusive measurement of the flattening of the Halo, we 
cannot place strong constraints on the merging rate onto the 
Milky Way. 

However, assuming that the M-giant stellar population 
in the Sagittarius stream is representative of other accre- 
tions, it appears that there has been no other more massive 
accretion than the Sgr dwarf {Msgr ~ 5 x 10* Mq) since 
that galaxy was accreted. The timescale for the accretion of 
the Sagittarius dwarf must be long, of order s everal Gyrs, as 
it takes of order 3-5 Gyr in the simulations (Ibata & Lewis 



199S; [bata et al. 2001a) for tidal tails to wrap around the 



Milky Way and appear as a band on the sky. 



5.2 Narrow streams from disrupted globular 
clusters 

As we show in Paper 1, the structures of particular interest 
are stellar streams from low mass systems such as globu- 
lar clusters. These give fine structures in phase-space which 
easily become dispersed if there is significant substructure 
in the Halo. No large scale stream other than that of the 
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Sgr dwarf was identified in the sample of late M-giant stars 
from the 2MASS 2IDR dataset above the 3a level. 

With our late M-giant sample, the 3a detection level 
corresponds to only 9 stars at the finest band resolution 
adopted (1°). The lack of detections reflect the fact that our 
M-giant selection criteria, which aimed to get rid of disk 
dwarfs, are very restrictive. A vast number of Halo K (and 
early M) giants were removed in this process leaving only a 
small number of Halo stars in the flnal sample. Applying our 
selection criteria to the 2MASS data of the massive (10^ M©) 
globular cluster Omega Centauri (distance 5.1 kpc) yields a 
total of 3 stars outside two core radii, or 5.8 stars correcting 
for the lack of data in the saturated central regions of the 
cluster. Only 1 in 2000 stars above the magnitude limit of 
K = 14.3 are kept after the colour-cut (though the stars 
that are retained are very luminous, and could be detected 
out to beyond 100 kpc). So it is not surprising that we have 
not found globular cluster streams in the first sample, eis the 
tidal stream of a completely disrupted globular cluster of the 
size of Omega Centauri would have insufficient late M-giant 
stars to bring the contrast of a polecounts peak above the 
3a detection level. 

The inclusion of Halo K-giants and early M-giants is 
needed to obtain high enough S/N to make a useful state- 
ment about dark matter dumpiness. Our second analysis, 
with a sample selected with practically no colour restric- 
tions, aimed to improve the statistics for narrow streams. 
However, no streams were detected above the 4a level, even 
though we had the capability of identifying narrow streams 
with as few as ~ 860 stars with the narrowest band we 
probed (0.5°). For comparison, Omega Centauri would pro- 
vide a stream of ~ 3450 stars in 2MASS, if it were located 
at a distance of 17 kpc, though the 25% sky coverage of the 
|6| > 30° region of the 2IDR dataset reduces this number to 
~ 860 stars. 

The detection of a dynamically cold thin stellar stream 
in the 2IDR would have called into question the existence 
of dark matter substructure. Our null detection is much less 
clear to interpret, however. We can rule out the existence 
of only very massive globular cluster streams (as massive as 
Omega Centauri), that are relatively nearby (< 17kpc), and 
that are confined to very narrow great circle bands (< 0.5°). 
This may mean that massive globular clusters are rare, or 
that the Galactic potential is significantly flattened at these 
radii, or that the potential has evolved substantially since 
the disruption of these stellar systems. 

Future experiments will help solve this question. By re- 
analysing the full 2MASS dataset, a ~ 3 times larger area of 
sky can be surveyed for streams (not ~ 4 times larger, since 
low latitude disk fields will likely be too complex to model). 
With that data in hand it may be possible to improve the 
stream detection limits by an order of magnitude, especially 
if the foreground Galactic component can be more accu- 
rately modelled in the absence of the complex 2IDR window 
function (which we believe was the primary obstacle in not 
reaching the Poisson noise limit in the polecount analysis). 
Other purely photometric experiments, such as the VISTA ^ 
survey or the PRIME ^ mission, will also greatly improve 



stream detection limits due to improved statistics resulting 
from a much greater photometric depth than 2MASS. Larger 
gains will come from the addition of astrometric information. 
A first step will be to use proper motions derived from pho- 
tographic plates (provided, for example by the GSC II 
with which it will become possible to obtain a cleaner Halo 
sample by requiring that distant Halo stars have proper mo- 
tions consistent with zero. The analysis of the datasets men- 
tioned until now can provide information on Halo streams 
only for the restricted case that the Galactic potential is 
nearly spherical, when streams are spatially-confined bands 
on the sky. The advent of GAIA (Ferryman et al. 2001), with 
its ability to accurately measure proper motions even for 
distant Halo stars, will permit the identification of streams 
even if the Galactic potential is substantially flattened, as 
the streams will be seen as phase-space clumps. 



6 CONCLUSIONS 

A pole-count analysis of Halo M-giants stars in the Sec- 
ond Incremental Data Release of the 2MASS dataset reveals 
a strong stream of stars along the orbit of the Sagittarius 
dwarf galaxy. The pole of the Sagittarius stream is located 
at exactly the po sition predicted by our previous analysis 
of Halo C-stars ([bata et al. 2001a), confirming the exis- 



tence of this large-scale stellar stream. This stream appears 
fairly narrow, ~ 12° in M-giants. The alignment of the nar- 
row stream along a great circle rules out the possibility that 
th e Halo mass distri bution is highly flattenned, as shown 
in [bata et al. (2001a ). Future radial velocity measurements 



of these M-stars will easily allow us to discriminate stream 
from Halo or foreground stars. Our conclusion regarding 
the shape of the dark matter Halo will be strengthed if the 
stream stars are found to be distributed all around the great 
circle track. 

Our main motive in undertaking this study was to iden- 
tify kinematically cold streams from low mass systems such 
as globular clusters, which could be used to trace the sub- 
structure of the dark matter Halo. However, no globular 
cluster streams were identifled in the sample of Halo late- 
type M-giants, though we argue that this is due to insuf- 
ficient statistics in the sample. Neither were any streams 
detected in the second sample that contained K and M- 
giants (as well as numerous foreground stars). In this sec- 
ond sample, nearby (< 17 kpc) globular cluster streams, of 
mass similar to that Omega Centauri, and that are closely 
confined to a great circle (< 0.5°) could have been detected. 
These stringent conditions greatly reduce the chances of de- 
tection of any Halo streams. In reality, even if the Halo is 
spherical, the flattened potential of the Galactic disk will 
cause some precession, especially to nearby streams. Also, 
the Sun's displacement away from the center of the Galac- 
tic potential means that only very distant streams can be 
located precisely on great circles. 

Better statistics are needed to allow the detection of 
globular cluster streams on wider great-circle bands, and 
deeper datasets are needed to probe more distant orbits. 
It will certainly be possible to improve this analysis in the 



sec it.t.p : / /hmu . vi ata . a r. . iil; 



see ittp : //prime .pha. jhu. edu/primemission. htm 



* see ^ittp : //www-gsss . stsci . edu/gsc/gsc2/GSC2home .htn 
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near future with the fuU sky 2MASS catalog, and better 
dwarf star rejection can be provided by proper motion se- 
lection (derived from the GSC-II catalog, for instance). In 
the longer term, the next generation of IR surveys (e.g., 
VISTA, PRIME), and the next generation of astrometric 
surveys (e.g., DIVA, GAIA) will greatly improve the num- 
ber of faint Halo sources and allow us to identify individual 
streams from their photometric and kinematic signatures. 



TuUy, B., Somervill e, R., Trentham, N., Verheijen, M. 2001, 
astro-ph/010753^ 
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Figure 1. The distribution of sources over the sky in the 2IDR dataset of the 2MASS project. Each 0.1° X 0.1° pixel (in rectangular 
coordinates) in the image displays logarithmically the number of sources in the corresponding region that have K < 14.3 (the ~ lOir 
limit of the photometry). 




Figure 2. The distribution in Galactic co ordinates of the 11 06 late M-giant Halo star candidates selected for analysis. The expected 
orbital track of the Sgr dwarf galaxy (from [bata et al. 2001a) has been overlaid. 
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I (degrees) 

Figure 3. The sky coverage for the present analysis. The black areas indicate the unobserved or unused regions of sky, while gray areas 
show the regions from which data were obtained, with a shading corresponding to the area subtended per pixel (0.01 square degrees on 
the Equator, decreasing to the poles). 
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Figure 4. The pole-counts diagram of our sample of Halo late-type M-giant stars, for a band of width 12.0°. The strong peak near 
{£ = 95°, b = 13°) reflects 15% of the sample, an ~ 5.5(7 overdensity. Its location, at the pole of the orbit of the Sagittarius dwarf 
galaxy, implies that we have detected the evolved stellar populations of the tidal stream of that galaxy, which was previously detected 
in C-st ars. T his polecounts diagram has a similar structure to the polecounts diagram constructed from Halo C-stars (Figure 5 of Ibata 
et al. 2C01e), with the difference that there is no evidence for the Magellanic stream in the present M-giant sample. This difference is 



likely a consequence of the lack of sky coverage around the South Ecliptic Pole in the 2IDR dataset. 



